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Abstract 



The Mathisson-Papapetrou-Dixon (MPD) equations for the motion of electrically neutral massive spinning parti- 
^ . cles are analysed, in the pole-dipole approximation, in an Einstein-Maxwell plane-wave background spacetime. By 
CN ; exploiting the high symmetry of such spacetimes these equations are reduced to a system of tractable ordinary differ- 
ential equations. Classes of exact solutions are given, corresponding to particular initial conditions for the directions 
7— ( . of the particle spin relative to the direction of the propagating background fields. For Einstein-Maxwell pulses a scat- 

^ ' tering cross section is defined that reduces in certain limits to those associated with the scattering of scalar and Dirac 

00 ' 

^ ' particles based on classical and quantum field theoretic techniques. The relative simplicity of the MPD approach and 
O ■ its use of macroscopic spin distributions suggests that it may have advantages in those astrophysical situations that 
^ involve strong classical gravitational and electromagnetic environments. 

(N : 
o 

^ 1 Introduction 

cr 

', Einstein's theory of gravitation predicts that certain non-stationary distributions of matter will produce dynamic 
^P' changes in the geometry of spacetime. Those that propagate are known as gravitational waves. Disturbances that 
J> . also involve accelerating sources of the electromagnetic field are also expected to produce electromagnetic waves. 
^ , Such scenarios can occur in many violent astrophysical processes and the energetics of such processes are poorly 
' understood in many cases. It has also been suggested that matter with spin may play a dual role in many astro- 
physical processes that emit electromagnetic and gravitational waves. The constituents of plasmas that participate in 
the generation of such waves may become polarised in strong magnetic fields and spinning matter may respond to 
tidal forces that include interactions of spacetime curvature with angular momentum. Such processes belong to the 
realm of general relativistic transport theory. However some insights into the interaction of spinning test matter with 
gravitational and electromagnetic waves can be gained by studying the classical motion of spinning test particles in 
plane gravitational and electromagnetic field solutions to the Einstein-Maxwell equations. In this paper we explore 
the motion of massive electrically neutral particles with spin in such fields. Since the electromagnetic field influences 
the spacetime metric it will have an effect on such electrically neutral particles even in the absence of gravitational 
waves. 

By neglecting self-gravitation and back-reaction the dynamics of classical test particles with angular momentum 
was first studied in detail by Mathisson, Fock, Papapetrou et al [^Jj. The theory was further clarified by Dixon [||, ^] 
using a rationalised multipole expansion technique and developed by Ehlers, Rudolph [f|] and others. In this article 
we employ the pole-dipole approximation so that the motion of a classical spinning test particles is governed by 
the Mathisson-Papapetrou-Dixon (MPD) equations. Such MPD equations predict gravitational spin-spin interactions 
between rotating stars and orbiting massive spinning particles. 

The properties of plane gravitational wave spacetimes have attracted considerable attention due to their high sym- 
metry [0]. This has led to investigations of colliding gravitational and electromagnetic waves [^] and the scattering 
of massive particles in such backgrounds The scattering of scalar fields was considered by Gibbons [§], Garriga 
and Verdaguer Jl0| ] and others. The corresponding differential scattering cross-sections were found to be similar to 
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those for the scattering of classical spinless particles. The scattering of electrically neutral Dirac fields was subse- 



quently investigated by Bini and Ferrari [ ]11| ] and compared with that involving classical and quantum spinless particle 
scattering. 

The classical scattering of neutral spinning particles by plane gravitational waves has been considered from a 



number of different perspectives [|12|, [T3| [Bp. In a recent paper QT3[ ] a special class of solutions of the MPD equations 
describing the non-geodesic motion of a massive spinning test particle in plane gravitational wave spacetime was 
constructed. The spin-curvature coupling was shown to give rise to parametric excitations of spinning matter by 
harmonic plane gravitational waves. 

In this paper we generalise this analysis and suggest that information about the scattering of electrically neutral 
spinning matter in an Einstein-Maxwell plane wave background can be gleaned from solutions to the MPD equations 
by exploiting the high Killing symmetry of such spacetimes. 

In section (2) the basic MPD equations are summarised and an expression given for constants of the motion 
that can be generated from Killing symmetries. Such expressions are pivotal in the subsequent analysis. In section 
(3) the MPD equations are expressed in different charts adapted to Einstein-Maxwell plane wave spacetimes. Such 
charts facilitate the the symmetry reduction of the equations to a tractable system of ordinary differential equations. 
In section (4) several exact solutions to these equations are discussed. The final section exploits the properties of 
an Einstein-Maxwell pulse (sandwich spacetime) in order to estimate the scattering of an arbitrary distribution of 
classical spins according to the MPD equations. The result is compared with the classical and quantum scattering 
cross-sections of scalar and spinning particles calculated by other techniques in such backgrounds. 

2 Equations of motion for neutral spinning test particles 

In the monopole-dipole approximation the MPD equations determine the world-line of a spinning test particle with 
tangent vector V, momentum vector P and spin 2-form s in terms of the spacetime metric g, its Levi-Civita connection 
V and Riemannian curvature tensor R. [] By introducing the metric duals p = P = g(P, — ) and v = V and denoting 
interior operator (contraction with any vector V) on forms by iy [16], the MPD equations ^ can be written in a 
compact form as 

V v p = i v f (1) 
Vys = 2p A v (2) 



with the spin condition [18] 



i P s = (3) 

and forcing term 

/ = ^(i^A^KAe 6 (4) 

in terms of the exterior product A and Hodge map * associated with a metric with the signature (—,+,+,+). Thus 
the "inverse Hodge map" denoted by [16] acts on any g-form oj according to -k^uj = — (— l)'?( 4 ~ l ?) satisfying 



*(* 1 uj) = * 1 {*uj) = u>. The curvature 2-forms i? a j, in any basis {X a } with co-basis {e a } are related to the 
curvature tensor components R a bcd by 2 R a b(X c , Xj) = R a bcd, (a, b,c,d = 0, 1, 2, 3). In terms of the components 
fa = f(X a ), v a = e a {V), p a = e a (P), s ab = s(X a ,X b ), equation (§ becomes 

r 1 d b cd 

Ja = ~2 RabcdV S . 



'Units in which c = G = 1 are used throughout. 
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It follows that [Q] the velocity of the spinning test particle is proportional to the tangent vector V whose components 
take the explicit form 

p c v c f a 2R bcde p c s ab s de \ 

' P ~ T-r: s ^7-77 • (5) 



pf p f \ 4p c p c - R bcde s bc s de 

The freedom to normalise V ensures that p c v c is arbitrary and once a parameterisation of the world-line has been 
chosen <^> permits the computation of the world-line PJ]. The norm of the momentum vector given by 



m = y/-g(P,P) (6) 
may be identified as the mass of the particle. This is used to define the normalised momentum vector given by 

P 

U =-. (7) 

m 



In Minkowski spacetime U coincides with the 4- velocity V/ \/—g(V, V) but in general these two vectors differ in 
curved spacetime. Using the 1-form u = U the spin 1-form I is defined by I = — | * (u A s) and the spin vector by 
L = I. The angular-momentum 2-form s may be related back to the spin 1-form by 

s = 2i u *- l l. (8) 



The norm of the spin vector L denoted I = yj g(L, L) defines the "spin" of the particle. The reduced spin vector and 
1-form are defined by S = L/m and ? = l/m respectively such that the norm a = g(E, S) = i/m denotes the 
spin to mass ratio. It follows from ([!]), (§) and (||) that a and m are constants of motion and the vectors U and S are 
orthogonal, i.e. g(U, S) = 0. 

A vector field K on spacetime is a Killing vector field if 

C K 9 = 
£ K F = 

for any Einstein-Maxwell solution {g, J 7 } satisfying 

G a + 8nT a = (9) 



d*J r = (10) 

where 

T a = i {% Xa T A^T -% Xa ^T AT) (11) 

in terms of the Einstein 3 -forms Q a in any basis {X a } [|l6|]. 

For each such K one may find a constant along the world-line C. Sufficient constants of the motion enable one to 
integrate the equations of motion in terms of properties of the background metric. Given K = k a X a it follows from 
([[]), @ and @ that the quantity 

C K = ★- 1 |xA*p+id^A*«| = k aP a -^k [a . b] s ab (12) 

is preserved along the world-line of a spinning test particle [^,15]. Its value can therefore be fixed in terms of 
particular values of p and s at any event on their trajectories. 
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3 The MPD equations in an Einstein-Maxwell plane-wave background 

The metric for an Einstein-Maxwell plane wave spacetime takes the form 

g = H{pL, x, y) d/t <E> d/t — — (d/2 ® du + dz> ® d/t) + dx ® dx + dy (g> dy (13) 

in "harmonic" (Kerr-Schild) coordinates (i,x,y,z). It describes a plane gravitational and electromagnetic wave 
travelling in the ^-direction j^, [7|], where fi = i — z, v = i + z and 

H(p, x, y) = g (A) (£ 2 - y 2 ) + em (A) (£ 2 + y 2 )- (14) 

Here 4> g (fi) is an arbitrary gravitational wave profile whereas </> em (A) is related to the electromagnetic potential 1-form 

A = Ai{fi) dx + A 2 {fi) dy (15) 
with arbitrary waveforms A\ (fi) and A2 (A) by 

<M£) = -2vr (^i(A) 2 + -4(A) 2 ) (16) 

and ' denotes differentiation. The electromagnetic field 2-form T = d A. 
In the orthonormal co-basis defined by 

e° = f + (l-#GM, 2 }))f 
e 1 = dx, e 2 = dy 

e 3 = ^-(l + tf(A,*,y))f (17) 

let {X a } be the dual basis satisfying e a (X b ) = 5 a b . This basis is parallel along the geodesic observer O : r 1— > 
(i(r) = r, x = 0,y = 0, i = 0), i.e. V X a \o = 0, and takes the form {Xq = d^X\ = 8x1X2 = d y ,Xs = dg} 
and sets up a "local Lorentz frame" along O. In these coordinates this observer curve is parameterised by O : r 1— > 
(/t( r ) = r, j>( T ) = r> £ = 0, y = 0). 

We are interested in the motion of a spinning test particle with mass m, spin £ = a m initially at "rest" relative 
to observer O and excited by a pulse of plane gravitational and electromagnetic waves . Thus we parameterise the 
world-line of this particle by C : A 1— > (/t(A) = A, z>(A), £(A), y(A)) for functions z)(A),x(A),y(A) subject to the 
initial conditions: 

i>(0) = 0, x(0) = x in , y(0) = y in (18) 
with arbitrary position parameters x- m and y m . The tangent vector V to the worldline C is written in these coordinates: 

V = C = d A + z)'(A)^ + x'(A)^ + y(X)dy. 

The pulse nature of the wave can be accommodated by taking the arbitrary profiles in the Einstein-Maxwell solution 
above to have compact support in the variable fi. Where such profiles vanish the spacetime is flat. 
Since the mass m and spin a of the test particles are constants of the motion one has 

Vab u a (X)u b (X) = -1 (19) 
Vab ^(X)^(X) = a 2 (20) 
Vab u a (X)s b (X) = (21) 
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in terms of the orthonormal components u a = e a (U) and ? a = e a (S) with rjab = diag(— 1, 1, 1, 1). The initial 
conditions for u a are: 

u°(0) = 1, u^O) = u 2 (0) = u 3 (0) = (22) 
and we write the initial conditions for s a 

,°(0) = 0, ^(0) = ,1, , 2 (0) = <&, , 3 (0) = & (23) 

where the constants ^ satisfy (<^) 2 + ) 2 + ) 2 = a 2 . 

For convenience let 

<M£) = 0g(£) ± 0em(£) (24) 
*±(A) = /%±(C)dC. (25) 

JO 

The Killing symmetry of the Einstein-Maxwell plane wave spacetime is most conveniently expressed in "group" 
(Rosen) coordinates (/j,,r/,x,y). For suitable ranges these are related to the previous chart by the relations: 

p, = fj, (26) 
f> = v + x 2 a(n) a'(fi) + y 2 b(ji) b\n) (27) 
x = a(fj,) x (28) 

y = b(fjL)y (29) 
where the metric functions a(n) and 6(/x) satisfy 

a"(p) = ^+(a*)o(a*) (30) 

= -</>-&) b(jj.) (31) 
with a convenient choice of the (gauge fixing) conditions: 

o(0) = 6(0) = 1 (32) 

a'(0) = 6'(0) = 0. (33) 
In this coordinate chart the metric g takes the form 

g = -- (d/i (g> du + dv <g) d^x) + a(/i) 2 dx <g> dx + b(fi) 2 dy ® dy. (34) 

This Einstein-Maxwell spacetime admits five independent Killing vectors : 

d d d 

Ki = —, K 2 = —, K 3 = — (35) 

ox ay 

K 4 = AGu) — + 2x— (36) 
ax ov 

X 5 = 5( M ) — + 2y— (37) 
ay ov 
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where 



2 ■ 



(38) 



In these coordinates we write the tangent vector 

V = d„ + z/(A)<9„ + x'(X)d x + y'(X)d y . 
It follows from ©, @, © and © that 



r??. 



Cjfj = — — , Cx 2 = 0, Ck 3 = 0, C^- 4 = m/3 2 , C^ 5 



-m/3i 



(39) 



where f3\ = <^ + y; n and = — ^in- From ( |T8| ) and ( |26| ) - ( |29^ ) in group coordinates the initial conditions for the 
world-line of the spinning test particle C : A \-* (/j(A) = A, u(X),x(X),y(X)) read 



i/(0) = 0, x(0) = Xi, 



2/(0) = 



(40) 



The five Killing vectors above yield five algebraic relations, (fl2|), between the components of u and q for any A. 
Equations (|l^), (p0|), ( |2l| ) provide three further algebraic relations. These eight equations may be solved in terms of 
the constants of motion and initial conditions to yield the solutions (f|l]) and (42) below: 



n (A) = 1 + ^ a ' (A) 2 + ^ 6 / (A) 2 



u\X) 



2 w 2 
foa'{\), u 2 (X)=p 1 b'(X) 



u*(\) = ^a'(Xy + ^b'(X) 



(41) 



and 



^°(A) 



{(/3 2 V(A) 2 + /3 2 &'(A) 2 ) , 3 (A) 



/3 2 a'(A) 2 + /3 2 &'(A) 2 -2 
- 2(y(A) - p 1 )p 2 b(X)a'(X) - 2(x(X) + p 2 )Pia(X)b'{X)} 

1 



= -(y(A)- fr)6(A) + 



{-2/? 2 a'(A) ? 3 (A) 



/3 2 a'(A) 2 + /? 2 6'(A) 2 -2 
+2 ft a'(A) [(y(X) - ft)ft6(A)a'(A) + (x(X) + ft)fta(A)6'(A)] } 



? 2 (A) = (x(X) + ft)a(A) + 



{2ft 6'(AX 3 (A) 



ft 2 a'(A) 2 + ft 2 6'(A) 2 -2 
-2 A 6'(A) [(y(A) - ft)ft6(A)a'(A) + (x(X) + ft)fta(A)6'(A)] } . 
Furthermore they imply a quadratic equation for ? 3 (A): 

? 3 (A) 2 + P(A) ? 3 (A) + Q(A) = 

where 

T{X) = -2 (y(X) - ft)ft6(A)a'(A) - 2 (x(A) + ft)fta(A)6'(A) 



(42) 



(43) 



(44) 
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and 



Q(A) = -^{ftV(A) 2 + ft 2 6'(A) 2 -2} 2 

+ 0/(A)-ft) 2 6(A) 2 | 2 ^ a / (A) 2^ (A) 2 + ^4 fl / (A) 4 + ^4^)4 _ 4^2 & / (A) 2 + 4 | 

+ (,(A)+ft) 2 a(A) 2 {20*0*^2^2 + ^ (A) 4 + ^ (A) 4 _ 4/3 2 a , (A)2 + 4} 
+2ftfta(A)6(A)a'(A)6'(A)(x(A) +ft)(y(A) - ft). 



(45) 



Thus for given metric functions a(A), 6(A) one has solutions ? 3 (A) in terms of x(X) and y(A). 

The differential equations for the components ^(A), x(A), y(A) of the tangent vector V are obtained from substi- 
tuting ©, (|), © and © into (|): 

x'(A) = 2 0_(A) {fta(A)6(A)6'(A)(x(A) + ft)(y(A) - ft) + ft6(A)V(A)(y(A) - ft) 2 

-(y(A) - ft)6(A) ? 3 (A)} / (a(A) (ft 2 a'(A) 2 + ft 2 6'(A) 2 - 2) } - (x(A) (46) 

y'(A) = 2 0+(A) {fta(A)6(A)a'(A)(x(A) + ft)(y(A) - ft) + fta(A) 2 6'(A)(x(A) + ft) 2 

-(.(A) + ft)a(A), 3 (A)} / {6(A) (ft 2 a'(A) 2 + ft 2 6'(A) 2 - 2) } - {yW '^'^ (47) 



v'{X) = 1 + (x(A) + ft) 2 a'(A) 2 + (y(A) - ft) 2 6'(A) 2 
<MA) 



+ 



+ 



ftV(A) 2 + ft 2 6'(A) 2 - 2 {" 2 °( A ) 2 + ^ (^«(A) 2 + 2ft &'(A) 2 y(A) - ft 2 6'(A) 2 - 2) 
4ft a(A)6(A)a'(A)6'(A) (x(X) + ft) (y(A) - ft) 2 + 4a(A)6'(A) ? 3 (A) (x(X) + ft) (y(A) - ft)} 
MA) {2 6(A) 2 (y(A) - ft) 2 (ft 2 6'(A) 2 - 2ft a'(A) 2 *(A) - ft 2 a'(A) 2 - 2) 



ft 2 a'(A) 2 + ft 2 6'(A) 2 -2 

-4ft a(A)6(A)a'(A)6'(A) (x(A) + ft) 2 (y(X) - ft) + 46(A)a'(A) ? 3 (A) (y(X) - ft) (x(X) + ft)} (48) 
These equations are decoupled as the sets: 



and 



(x(X)\ 

\y(x)J 



F 1 (x(X),y(X),, 3 (X)) 
F 2 (x(X),y(X),^(X)) 



(49) 



where from (4-3), 



u'(X)=F 3 (x(X),y(X),, 3 (X)) 



, 3 (X) = F 4 (x(X),y(X)). 



(50) 



Thus the world-line of the spinning particle can be determined by solving the two coupled non-linear differential 
equations (49) for x(X) and y(A) and inserting these solutions into (|50|) and solving for v{X). The evolution of the 
reduced spin vector £ then follows by solving (42) and (43). 

Given the initial position parameters x- m , y m and reduced spin vector components <q^, <;?, the evolution of the 
normalised momentum vector U follows more simply from ( f4"l| ) for arbitrary profiles 4>±(fi) and the corresponding 
a(ji), b{n) obtained from ©, ©, |||), ©. 



4 Particular solutions describing geodesic and non-geodesic motions 

Since the above equations are non-linear it is difficult to find general analytic solutions for the particle's worldline 
and spin. It is however possible to find interesting particular analytic solutions for special initial conditions. Thus 
if the initial spin vector is in the direction of propagation of the gravitational and/or electromagnetic wave then it 
remains so and furthermore the subsequent motion is geodesic. With ^ = = 0, a = |q^| and any x m and y m , the 
world-line C is simply: 

/x(A) = A, x{X) = xin, y(X) = y in , z/(A) = A (51) 
with normalised momentum components: 

u°(X) = l + ^a'(X) 2 + yfb'(X) 2 
u^X) = xi„ a' (A), u 2 (X)=y in b'(X) 

u 3 (X) = ^ a '{X) 2 + yfb'{Xf (52) 

and reduced spin components: 

? (A) = -|(4a'(A) 2 +^&'(A) 2 ) 

^(A) = - ? f n x in a'(A), e{X) = -&y m b'(X) 

^(X) = 'f(l-x 2 n a'(X) 2 -y 2 n b'(X) 2 ). (53) 

A class of solutions describing non-geodesic motions can also be found in which the particle spin vector is 
transported V-parallel along its world-line, transverse to the direction of wave propagation. Consider the case with 
a = |«tal.«ta = <fa = Oand 

x(A) = x- m . (54) 
c 3 (A) = x in + y in - y(X)) b(X)a'(X) (55) 



Then ( |46| ) implies 
and ([f7|) becomes 

with the solution 



m = (Aiv^vm (56) 



y(A) =ci(l-^y) + '/■„ (57) 



8 



satisfying the initial condition y(0) = y m . Equation (48) then reduces to 

z/(A) = l + (4) 2 fo-W + ^f) 
Using ( f3"ij ) and the initial condition u(0) = this can be integrated : 

y(A)\ 



KA) = A+(^ *-(A) 



6(A) ; • 



Substituting the above equations into ( |41| ) and ( |42| ) gives 



2 / 1 \ 2 

n o (A) = 1 + ^ a '(A) 2 + ^ + yia > b'(X) 2 
u\X) = x in a'(A), u 2 (X) = + y in ) 6'(A) 



« 3 (A) 



(4 + Vi 



'"' t'(A) 2 



and 



? °(A) = ^x in a'(A), ? 1 (A) = ^, ? 2 (A) = 0, ? 3 (A) = d^a'(A). 
Using (26) - ( |29| ) solutions (|54|), (|7]), ( |59| ) in harmonic coordinates become 



XT, 



(4 + Kh 



f(A) = A + ^ o(A) a'(X) + ^y*; 6(A) 6'(A) - 2 ^ + y in ) 6'(A) + ' 
x(X) = a(X) x in , y(A) = 6(A)(? i 1 n + y in ) - <^ n 



1 \2 



$_(A) 



(4 + 2/in)' 



1\2 



z(A) = -^o(A)a'(A)+ , 



6(A) 6'(A) - 2 + y m ) 6'(A) + M- *_(A). 



(58) 



(59) 



(60) 



(61) 



(62) 



For such non-geodesic motion the difference between the corresponding momentum and velocity vectors is nonzero 
and takes the form: 



U-V 



2(4r<MA) 



d_ 

dv 



(63) 



With the alternative initial spin conditions = = and a = |^|, one similarly finds 

n o (A) = 1 + (^-^) 2 q/(a)2 + vl_ 6/(A)2) 



u x (A) = -( ? 2 -x in )a'(A), u 2 (A) = y in 6'(A) 



(64) 



and 



?°(A) = 4yin&'(A), ? X (A) = 0, ^ 2 (A) = ? 2 , ? 3 (A) = ? 2 y m 6'(A). 



(65) 
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The corresponding world-line is then given by 

f( A ) = A ^""^ 2 a(A) a'(X) + V f 6(A) b'(X) -2&(& - x in ) a'(A) - ^ * + (A) (66) 

x(A) = -a(A)(4 - x in ) + q? (67) 
y(A) = 6(A)y in (68) 

5(A) = ~ Xi " )2 a(A) a' (A) + V f 6(A) 6'(A) - 2 ^ - x to ) a'(A) - ^ $+(A) (69) 
in harmonic coordinates. In this case the difference between momentum and velocity vectors is: 

U-V = 2(4)V f (A)JL (70) 

5 Scattering cross-section due to an Einstein-Maxwell wave pulse 

Aside from regularity conditions the above calculations have not imposed any strong conditions on the structure 
of the profile functions that characterise the Einstein-Maxwell background spacetime. Consider now the scattering 
of neutral massive spinning particles by a plane gravitational and electromagnetic wave pulse represented by the 
metric (|l3|) with functions 4> g (fi) and 4> em (£i) that vanish for fi < and fi > 5, where 5 denotes the width of the 
wave pulse. Such a sandwich wave gives rise to two half Minkowski spacetimes namely M.[ B for fx < and Ai out 
for ft > 5 corresponding to the domains before and after the passage of the plane wave pulse. Within each half 
Minkowski spacetime the orthonormal basis {X = d^, X\ = d&, X 2 = dy, X% = dg} defines a Lorentz frame, and 
{X a } along the geodesic observer O still remains a local Lorentz basis also in the "sandwiched" domain V where 
< fi = i — z < 5. Consider a collection of spinning particles, labelled by their initial x and y components, i.e. 
x m and y- m . Suppose these particles are initially at rest relative to O having common initial velocity and momentum 
given by V(0) = U(0) = Xq but with possibly non-uniform initial spin orientations described by: 

?in = finO^iflj 2/in)j ?in = fin O^in ; Um)i Sin = ^in^in) 2/in) • (71) 

Particles scattered by the pulse will have different out-going momentum vectors, depending on their initial data. In 
A4 m and A4 out these vectors have different constant components u a and changes in these quantities along C across 
the domain V reflect their gain in energy-momentum in the Lorentz frame associated with O. In any chart for Ai m 
write a' out = a' (A), 6o Ut = 6'(A) and p% ut = m u a (X) for any A > 5. In such a sandwich spacetime the components of 
the final momentum vector are readily evaluated from ([fj): 

Pout = ™ + j(&- x- m ) 2 a' oat 2 + j (<£ + y in ) 2 b' out 2 (72) 

Pout = ~ m (<5n - x in) a'out ( 73 ) 
Pout = m (4i + ^n) b 'out (74) 
Pout = J - ^in) 2 Gout 2 + -J + Mn) 2 b 'oui ■ ( 7 5) 



Equations ( |73| ) and ( |74[ ) define a map between the transverse components of momentum in any spacelike plane 
associated with O's rest space in -M ut and the initial location of particles in a similar plane in A4- m . It is therefor 
natural to define a classical differential scattering cross-section associated with this map in terms of the out-going 
momentum components (pl ut ,Pout)'- 

, _ , , , a'PoutflPout f, , dsl d<& d<^ d& \ _1 

da = dx in ay in = — =— - — tj— 1 + ^ j h -5 — -5 j. — ~f. — • (76) 

m a outOout V <%n ox in dx in dy m dy in dx in J 
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This cross-section provides a means of estimating the effects of a non-uniform distribution of spin on the scattering 
of particles by a passing Einstein-Maxwell pulse of radiation. When this initial polarisation is uniform (i.e. are 
independent of the parameters x- m , y- m ) this expression agrees with the differential cross-section describing the scatter- 
ing of classical or quantum non-spinning particles JTo| ] by an Einstein-Maxwell pulse and furthermore coincides with 
the differential cross-section describing the scattering of polarised incident quantum Dirac particles in such a back- 
ground [pH]. A priori there is no reason to expect that electrically neutral particles with spin-curvature interactions 



described by the MPD equations should behave in this way. More generally we believe that the cross section (76) de- 
scribes how such particles with a non-uniform spatial distribution of spins are scattered by an Einstein-Maxwell pulse 
in space. In the absence of a relativistic transport description of astrophysical phenomena involving spin-curvature 
interactions the approach above offers a relatively straightforward mechanism to estimate the significance of such 
interactions. 



6 Conclusions 

By exploiting the high symmetry of Einstein-Maxwell plane-wave spacetimes we have reduced the pole-dipole MPD 
equations for the motion of massive particles with spin to a system of tractable ordinary differential equations. Classes 
of exact solutions have been found corresponding to particular initial conditions for the directions of the particle 
spin relative to the direction of the propagating background fields. For Einstein-Maxwell pulses we have defined 
a scattering cross section that reduces in certain limits to those associated with the scattering of scalar and Dirac 
particles based on classical and quantum field theoretic techniques. Such techniques are considerably more intricate 
than those discussed here. The relative simplicity of our approach and its use of macroscopic spin distributions 
suggests that it may have advantages in those astrophysical situations that involve strong classical gravitational and 
electromagnetic environments. 
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